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Abstract 

This work proposes a novel technique for the numerical calculus of 
" ^ the fractal dimension of fractal objects which can be represented as a 

^> closed contour. The proposed method maps the fractal contour onto 

i-C a complex signal and calculates its fractal dimension using the Fourier 

transform. The Fourier power spectrum is obtained and an exponen- 
tial relation is verified between the power and the frequency. From 
t— I the parameter (exponent) of the relation, it is obtained the fractal 

dimension. The method is compared to other classical fractal dimen- 
sion estimation methods in the literature, e. g., Bouligand-Minkowski, 
box-couting and classical Fourier. The comparison is achieved by the 
calculus of the fractal dimension of fractal contours whose dimensions 
are well-known analytically. The results showed the high precision and 
robustness of the proposed technique. 

(N 

1 Introduction 

In the last years, the fractal dimension has become an important tool for 
characterize objects, forms and surfaces in some areas of the science. The 
examples of fractal dimension applications are in a wide range and can be 
found in areas so distinct as medicine |Tj, texture analysis [2J, geology j3], 
botany [HE], materials engineering [6], electronic [7|, physics [8], histology 
[9] , soil analysis [10] , plant diseases [11] , polymer analysis [12] , among many 
others. 
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The importance of the fractal dimension and consequently its large num- 
ber of applications can be explained by the fact that this measure gauges, 
at different scales, how the fractal fills the space in which it is immersed. 
Thus, the fractal dimension characterizes the fractal object with a number 
which is not dependent from the observation scale. As a consequence of 
this multiscale characteristic, this measure can be still used to detect the 
level of self-similarity in the analyzed fractal object, that is, it shows how 
much the object is composed by similar copies of itself at different scales, an 
intrinsic characteristic of fractals. The self-similarity property is also widely 
observable in natural objects [13J as the leaves of a tree or the stream of a 
river and this fact strongly encourages the researchers in using the fractal 
dimension to characterize objects from the real world, as in the applications 
previously cited. 

The fractal dimension concept was originally defined in |13j as being the 
Hausdorff dimension [ll] of the fractal object. This concept of dimension is 
an analytical measure whose calculus uses to be complex and even impos- 
sible in many cases. Indeed, as the Hausdorff dimension is calculated by a 
covering process using infinitesimal sets, it is necessary for this calculus to 
know a well defined mathematical rule used in the formation of the fractal. 
Such need is not satisfied for example in some Lindenmayer fractals as the 
three-dimensional weeds [15J. As a consequence, along the years, a lot of 
numerical methods were developed in order to estimate the fractal dimen- 
sion in these cases. These methods became still more important with the 
representation of fractal objects in digital images. In the digital representa- 
tion, the fractals are shown in a discrete space (image grid) which derails the 
Hausdorff dimension calculus once the fractal is handled as a disconnected 
set of points, compromising seriously the resultant dimension. 

The approximated fractal dimension estimation methods can be split 
into two approaches. The first are the spatial methods, as box-counting 
[16], Bouligand-Minkowski [IT], mass-radius [18], dividers [19J, among oth- 
ers, based on geometric strategies used for measure the self-similarity of the 
object. The second approach is based on spectral measures, as Fourier [20] 
and wavelets f2~T] . Particularly, fractal dimension estimative methods based 
on Fourier are interesting once these methods have easy implementation and 
low computational cost using Fast Fourier Transform (FFT). Another ad- 
vantage of such techniques is the fact that they are less sensible to geometric 
transformations such as scale, rotation and translation [20j[22 . 

This work proposes a novel method based on Fourier fractal dimension 
to calculate the fractal dimension of contours. In this technique, the fractal 
object must be representable by a closed contour and this contour is mapped 
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onto a complex signal. The power law of the power spectrum is used to es- 
timate the fractal dimension using a strategy similar to that described in 
|23j . The proposed method is compared in the calculus of the fractal dimen- 
sion with other fractal dimension methods. Although the work presented 
here is a theoretical one, and just mathematical fractals was considered, 
the method can be used to estimate the fractal dimension of signals (in the 
real or complex domain) and in this way, can be used in a wide range of 
applications. 

This paper is divided into 7 sections, including this introduction. The 
following discusses the theory of fractal dimension. The third addresses the 
Fourier fractal dimension. The forth describes the proposed method. The 
fifth explains the experiments. The sixth shows the results of the experi- 
ments and the last section does the conclusion. 

2 Fractal Dimension 

Fractal objects were studied by mathematicians like Cantor and Peano, since 
the 19 th century, when they were known as monster curves. Such strange 
structures did not obey the rules of the traditional Euclidean geometry. 
For instance, in 1904, Helge Von Koch presented a curve composed by an 
infinite sequence of equalizer triangles and which showed finite area and infi- 
nite perimeter. Such condition represented a contradiction to the Euclidean 
classical rules [24"] . 

Although these objects brought a certain interest in the principle, it was 
only with the arising of computers in 1970 decade that the fractal theory 
assumed a fundamental role. In that years, Mandelbrot |13j was the first 
author to formalize the concept of fractal. He defined fractals as organized 
structures which have infinite complexity, present self-similarity (in some 
level) and whose Hausdorff dimension exceeds strictly the topological (Eu- 
clidean) dimension. 

The self-similarity suggests that a worthy measure in the characterization 
of fractals must present a unique value independent from the observation 
scale of the object, once we observe the geometric aspect preserved when we 
analyze the fractal more and more microscopically. The measure satisfying 
this requirement is the fractal dimension. The Hausdorff dimension [T5j used 
in the definition of Mandelbrot is an example of such measure which captures 
the whole range of scales by covering (feeling) the fractal with infinitesimal 
sets. 

For the definition of the Hausdorff dimension, we must initially define 
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the s-dimensional Hausdorff measure H S (F) of a subset F G K": 

H S (F) = \im H°(F), 

£-!>0 

where 

H°(F) = inf |^ \Ui\ s : U { is an e-cover of f| , 

in which || denotes the diameter in K™, that is, \U\ = sup\x — y\ : x,y £ U. 
The Hausdorff dimension properly is obtained through the following expres- 
sion: 

dim H (F) = {s}\ inf {s : H S (F) = 0} = sup {H*(F) = 00} . 

A more practical definition of the Hausdorff fractal dimension can be 
provided from a generalization of the concept of topological dimension ap- 
plied to Euclidean objects. 

In fact, in Euclidean geometry, if we have an object A of dimension 
D and an object B which is the same object A, but with its linear size 
reduced in all the D directions by 1/7, the number of objects like B needed 
to cover (fill) the object A is given by N(l) = I . Notice that this formula 
is not dependable from the absolute linear size of A and B. The Figure [T] 
illustrates this fact in Euclidean geometry. In this way, the dimension may 
be calculated by: 

D = loam)) 
log{\) 

Extending to fractal structures, we have the following simplified version of 
the Hausdorff dimension: 

D = lintog (iyw) i 

«-° loud) 

in which N(e) is the number of objects with linear size e needed to cover the 
whole object |14j . For example, in the Koch curve, we have at each iteration, 
the curve covered by 4 segments whose length is 1/3 of that segment in the 
anterior step. The Figure [2] illustrates the process. Applying the dimension 
formula for each iteration k, we have: 

D = Ho, = Urn = „ L26 . 

6->0 l og {\) k^log(^) k log(3) 
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Figure 1: The relation between the scale factor L and the number N of 
objects covering the original structure in the dimension D. 



Figure 2: The first iteration in the composition of the Koch snowflake curve. 
4 segments with one third of the length of the anterior segment replacing it. 

The Hausdorff dimension in the most of cases presents a complex math- 
ematical calculus and in many cases this calculus is impossible, like when 
we do not know the rules of composition of the fractal. In order to fill this 
gap, along the years, a lot of numerical methods have been developed for the 
calculus of the fractal dimension in general situations. These methods can 
be split into two approaches: the spatial and the spectral techniques. The 
first is based on spatial measures as box-counting |16j , Bouligand-Minkowski 
[T7] . mass-radius |18j , dividers [19j . among others. The second approach is 
based on operations on the frequency domain, as Fourier |20j . wavelets |21j . 
etc. 

From the spatial methods, we explain briefly the operations of the two 
most used in the literature, that is, box-counting and Bouligand-Minkowski. 

2.1 Box-counting 

For the definition of the fractal dimensions in 3? 2 here presented, we consider 
the fractal object as a nonempty subset in 3ft 2 . The essential definition of 
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the box-counting dimension dimsC of C is therefore given by: 

dimsC = lim — ^- — -, 
<5->o —logo 

where N$(C) is the minimum number of sets with diameters at least 5 needed 
to cover C. With the aim of simplifying the numerical calculus, we can 
consider Ng(C) as the number of (5-mesh boxes b which intersect C, where 
b is a 5-coordinate of 3ft 2 given by: 

b=[b 1 5, (61 + 1)5 x b 2 5, (b 2 + 1)6], 

where b± and 62 are integers. 

For practical implementations, the box-counting method is based on the 
division of the space of the object to be analyzed into a grid of squares 
(boxes). At each iteration of the method it is used a different side length r 
for the boxes in the grid. The fractal dimension is given by —a where a is 
the slope of a straight line fitted to the curve of the function N(r) x r in a 
log-log scale, where N(r) is the number of boxes of length r which intersect 
the fractal object. The Figure [3] illustrates the process. 

2.2 Bouligand-Minkowski 

In the case of the Bouligand-Minkowski method, we initially must define the 
Bouligand-Minkowski measure of C given by: 

V(dX g rS) 
meas M (C,S,T) = hm — , 

r->-0 r 

where S is a structuring symmetrical element with radius r and V is the 
volume of the dilation between S and the boundary DC of C. The Bouligand- 
Minkowski dimension depends only upon the structuring element used and 
it is thus obtained through: 

dimM(C, S) = inf {r, measM{C, S, r) = 0} . 

For practical purposes, the Bouligand-Minkowski dimension is calculated 
through neighborhood techniques. Therefore, in cases where we analyze 
curves in the 5i 2 as is the case in this work, each point of the curve is 
replaced by a disk S e with diameter e composing a dilated area for each 
value of e. The fractal dimension in this process can be calculated by the 
following expression: 

^ M( c) = iim ( 2 _ iQ9Mxes e 

e->o y loge 
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Figure 3: The grid used in box-counting method with different size of each 
square and the squares intersected by the contour shaded in each case. The 
square size is represented by r and the number of intersected squares by N. 
(a) Grid with (r = 4). (c) Grid with r = 2. (d) Grid with r = 1. (d) Table 
of r by TV and graphical of log(r) x log(N) fitted by a stright line and its 
slope. 
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Figure 4: Bouligand-Minkowski method. A curve surrounded by circles 
with different radii composing the dilation area (shaded). The dilation is 
composed by the union of the circles and inasmuch the radius is increased, 
the circles begin to interfere in each other. The interference pattern follows 
the multiscale properties of the fractal allowing the calculus of the fractal 
dimension. 

In the computational calculus, the Bouligand-Minkowski method is ap- 
plied by constructing circles with radius r centered at each point of the 
fractal. The radius varies at each step of the method and the set of points 
pertaining to the union of the circles of radius r composes the dilation area 
A(r). With the aim of optimizing, the dilation area is calculated usually 
by a method named Euclidean Distance Transform (EDT) [5S].The fractal 
dimension is given by N — a where a is the slope of a straight line fitted to 
the curve of the function A{r) x r in a log-log scale and N is the topological 
dimension of the Euclidean space in which the fractal object is immersed, 
in this case, 2. The Figure [4] illustrates the method. 

Given a brief explains of the most commonly used fractal dimension 
spatial methods, this work will focus on the frequency approach and the 
next section describes the Fourier fractal dimension calculus. 

3 Fourier Transform and Fractal Dimension 

The use of the Fourier transform as a tool for the calculus of the fractal 
dimension is related to the definition of fractal transform. For our study, we 
must consider the Fourier transform of a mass distribution fi(u), u S K n , 
that is, a measure defined on a bounded subset of 3? n , for which < m(3? n ) < 
infty. The Fourier transform can thus be obtained through: 



where • denotes the scalar product defined over 9? n . 

From the Newtonian physics, we obtain the gravitational s-potential p s , 
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s > 0, of the mass distribution by: 

Ps(x) = I , — —^dm(y) 
J \x-y\ 

Extending the Physics analogy, the s-energy e s is obtained through: 

e s (m) = J p s (x)dm(x). 
By applying the Fourier transform to p s and using the Parseval's theorem: 

e s (m) = (2n) n c f 1(p s )(u)1(m(u))du, 
where x is the complex conjugate of x, so that: 

e s (m) = (2vr) n c J \u\ s ~ n \m{u)\ 2 du. (1) 

If there is a mass distribution m(u) on the set C £ S 2 for which the 
expression [l] is finite for some value(s) of s, so the Hausdorff dimension of C 
has its lower limit in s. Particularly, if |T(m(u))| < 6|u|~*/ 2 , for a constant 
value b, then e s (m) always converges if s < t. The greatest t for which there 
is a mass distribution m on C is called the Fourier fractal dimension of C. 

Russ |20j extended the use of Fourier transform to calculate the fractal 
dimension of objects represented as surfaces in a digital image. Russ associ- 
ated the Newtonian potential to the power spectrum of the discrete Fourier 
transform of the image. Based on this assumption, [20j[23j[26| shows that 
there is a simple exponential relation between the power spectrum of the 
Fourier transform of a gray-level image and the frequency variable in the 
same transform. Such fact can be expressed by the following formula: 

P(f) oc f-P 

It is still verified that the fractal dimension of the surface represented in that 
image is directly related to the coefficient (5 by the following expression: 

/3 = 2F + 2, 

where H is the Hurst coefficient. This coefficient is directed related to the 
fractal dimension FD by: 

FD = 3 - H 
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In this way, one may easily deduce that the fractal dimension of the surface 
can be obtained through: 



In practice, the Fast Fourier Transform (FFT) algorithm is applied to 
the image analyzed. The average power spectrum is taken along the fre- 
quency (radially) and the curve of the power spectrum against the frequency 
(distance from the center of the spectrum) is plotted in a log-log scale. A 
straight line is fitted to the curve and its slope is taken as the value of /3 
used in the formula above. 

It is also important to notice that for the best efficiency of this method 
in the calculus of the fractal dimension, the phase at each point of the 
transform must be the most random possible |2Uj . This fact implies that 
the values of the phase must be totally de-correlated, condition that may be 
checked by the calculus of a correlation measure along the Fourier phase. 



4 Proposed Method 

This work proposes a numerical method for the calculus of the fractal di- 
mension of fractal objects which may be represented by a closed contour. 
The contour is thus represented by a parameterized curve. Such strategy is 
largely used since it eliminates redundancies intrinsic to the binary image 
representation of the contour and identifies more explicitly the most impor- 
tant elements. The general approach to extract the parametric curve is to 
choose an arbitrary initial point pertaining to the contour and to walk along 
this contour following a pre-defined direction. Each point visited in the step 
t is the coordinate (x(t),y(t)) of the parametric curve. More details and 
methods for extract the parametric function of a contour may be seen in 

Posteriorly, following the order of the obtained contour, each point co- 
ordinate (x(t),y(t)) is mapped onto a complex-valued signal C{t): 

C(t)=x(t) + iy(t), 

where t is the parametric variable of the contour and i is the imaginary 
number. The Figure [5] illustrates the process. 

Thus, the Fourier Transform is applied to the signal C(t), obtaining the 
transformed signal C(u): 

/+oo 
C(t)e- 2nitu dt, 
-oo 
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Figure 5: A contour and its parametric representation, a) The contour of 
the snowflake fractal, d) Signals obtained from the x and y coordinates of 
the parametric curve. 



where u is the frequency variable. As for the classical Fourier dimension 
described in the Section |3j the function C(t) can be handled as a mass 
distribution and like in [20] the energy is obtained through the spectrum P: 

P = \%{C(u))\ 2 , 

where |.| denotes the conventional norm in K . As a consequence, the power 
spectrum becomes a powerful tool for the calculus of the fractal dimension 
of the object represented by the signal. The experiments verified that the 
generalization of Russ also holds in this case and a linear relation is observed 
between P and u, that is: 

Pocu D , 

where D is the fractal dimension. From the above expression, we can obtain 
the following relation for the dimension: 

D oc log u P. 

Applying the well-known basis change property and replacing the propor- 
tionality by using constants, we obtain: 

log(u) 

where k± and ki are real constants which do not affect the exponential 
relation. For a precise numerical calculus of D we must make use of intervals 



11 



with the smallest possible length and calculate: 

log(P + h)-log(P) 

D = ki hm - — - — — — — + k 2 . 

h->o log(u + n) — log{u) 

Disregarded the constants, the above expression represents the derivative of 
log(P) relative to log(u). And so we deduce that the fractal dimension can be 
simply obtained from the calculus of this derivative. For the calculus over the 
discrete signal C(t) the derivative is calculated in the practice by the slope 
of a straight line fitted to the curve of the graphical of log(P(u)) x log(u). 
From the experiments, we obtained the following values for ki and ki: 



Am 



k> 2 



1 

4" 



We also notice in the experiments a more significant variation of this 
slope depending on the region taken into account, as one can see in the 
Figure |6| 




log(frequency) 



Figure 6: The curve of power spectrum x frequency and the slope varying 
depending on the region of the graph to be considered. 

In the particular case of the fractal objects analyzed in this work, the 
experiments demonstrated that the slope has a great variation in the first 
points of the curve and after a certain number of points, the slope tends to 
remain constant. The Figure [7] illustrates such behavior by depicting the 
slope as a function of the first point taken into consideration for the linear 
fitting. 

Such fact was previsible since the fractal objects contain significant infor- 
mation in its more microscopic (detailed) scales. This interesting property 
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Figure 7: A graph representing the slope of the log-log power spectrum curve 
in function of the starting point considered in the graph. 

isnot observable in an Euclidean object but is intrinsic to fractals and is 
justifiable by the complexity inherent to fractals. The Fourier transform 
preserves the global aspect of the object in the first elements and the details 
in the last elements. The Figure [8] illustrates this fact by showing a fractal 
object (a), its reconstruction using the first elements of the Fourier trans- 
form (b) and the using last elements (c). We observe that the first elements 
preserved the global shape of the fractal while the last ones highlighted the 
details as the recursive entrances in the fractal shape. In this way, it is 
convenient the discard of the firt terms in the Fourier transform. 

For the application developed here, the first 10% of the points in the 
transformed curve were discarded, with the aim of capturing the maximum 
possible of details present in a fractal object. Such fact turns the method 
here exposed very sensible to imperfections in the structure of the fractal 
measured. 

5 Experiments 

For the validation of the proposed technique, this work calculates the fractal 
dimension of fractal objects whose Hausdorff dimension is well-known in the 
literature. 

The fractals used in the experiments are examples of fractal closed 
curves, which can be represented by a parametric contour. The curves used 
are Dragon curve [28] (Figure [9^a)), Fibonacci curve [13] (Figure [9^b)), 
Peano-Gosper curve [T3] (Figure [9^c)), the boundary of three variants of 
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(a) (b) (c) 

Figure 8: Fractal object reconstructed after Fourier transform, a) Original 
fractal, b) Reconstruction using the first terms, c) Reconstruction using the 
last terms. 



Julia sets (Figure |9^d) , (e) and (f)), Koch snownake curve [30] (Figure 
[g)), Douady- Rabbit curve |~9] (Figure ^h)), Terdra n curve [28] (Figure 
1)) and Vicsek curve [31] (Figure j9^j)). 
The curves had its fractal dimension calculated by the method proposed 
here and the result was compared to the analytical Hausdorff dimension 
shown in the literature and to the fractal dimension calculated by other 
numerical methods well known, that is, Bouligand-Minkowski, box-counting 
and Fourier spectral method. 

For a more robust test, the fractal curves were still submitted to rotation 
(by 45 and 90 degrees, representing a diagonal and a vertical perspective of 
the object), translation (by 10 and 20 pixels) and scale (by 0.5 and 2 times, 
representing a reduced and an extended copy of the original object). In 
a second moment, each fractal object was affected by a random punctual 
noise, as illustrated in the Figure 10 and the fractal dimension was measured 
by the compared methods. Such random noises are commonly detected in 
objects represented in digital images due, for example, to a scale reduction. 
It is difficult to remove them from fractal objects once the noise reduction 
methods use to smooth the curve, reducing in this way the details which are 
fundamental in the representation of fractals. Thus, for practical purposes, 
it is important for the fractal dimension calculus method to be robust to 
such noise. 

The accuracy of each fractal dimension method is measured by the aver- 
age modular error of each calculated value relative to the analytical Haus- 
dorff dimension value. 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 




Figure 9: Fractal objects used in the experiments, a) Dragon curve, b) 
Fibonacci curve, c) Gosper curve, d) Julia set 1. e) Julia set 2. f) Julia set 
3. g) Koch snownake curve, h) Rabbit curve, i) Terdragon curve, j) Vicsek 
curve. 




Figure 10: To the left, an example of a fractal and to the right, the same 
fractal affected by the punctual random noise used in the experiments. 
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6 Results 



The following tables present the fractal dimension calculated for the fractals 
affected by geometric transformations and noise. We can notice that the 
proposed technique is the unique totally invariant to the geometric trans- 
formations employed. Such fact is due to the intrinsic geometric invariance 
of the Fourier transform and to the manipulation of real values, avoiding 
rounding errors. 

The Table [T] shows the fractal dimension calculated by each compared 
method and the Hausdorff dimension. The tables [2] and [3] show, respectively, 
the fractal dimensions of the fractal objects rotated by 45 and 90 degrees. 
The tables [4] and [5] show, respectively, the fractal dimension of the fractal 
objects scaled by a 0.5 and 2 factor. The tables [6] and [7] show the fractal 
dimension of the fractal objects translated, respectively, by 10 and 20 pixels. 

The Table[8]shows the average error (mean of the differences between the 
calculated and the waited value) of the fractal dimensions calculated by each 
compared method. The Table [9] present the fractal dimension calculated for 
the fractals affected by the punctual noise. The Table 10 shows the average 
error of the fractal dimensions calculated by the methods. 



Table 1: Fractal dimension of original fractals for each compared method. 
The underlined values are that which presented the best approximation to 
the theoretical Hausdorff dimension. 



Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdorff 


Dragon 


1.4414 


1.4524 


1.3689 


1.5760 


1.5236 


Fibonacci 


1.3044 


1.3453 


1.2669 


1.2499 


1.2465 


Gosper 


1.1284 


1.1217 


1.0685 


1.1213 


1.1292 


Julial 


1.0905 


1.0538 


1.1283 


1.1485 


1.0812 


Julia2 


1.1465 


1.1972 


1.2146 


1.2509 


1.2683 


Julia 3 


0.9677 


1.0166 


1.0680 


1.0076 


1.0043 


Koch 


1.2479 


1.2497 


1.1080 


1.2722 


1.2619 


Rabbit 


1.2724 


1.3643 


1.2197 


1.3586 


1.3934 


Terdragon 


1.1871 


1.0949 


1.1162 


1.2313 


1.2619 


Vicsek 


1.3757 


1.4356 


1.3308 


1.5436 


1.4650 



From the Table [T] we observe that the proposed method was the most 
precise. It approximated more faithfully the dimension of the most of the 
fractals analyzed. In its turn, the experiments involving classical geometric 
operations (rotation, translation and scale), depicted in tables from [2] to [7J 
showed the robustness of the method. In the experiments, the fractal di- 
mension calculated by the proposed technique was completely invariant to 
these operations. This invariance is easily explained by the intrinsic geo- 
metric invariance of the Fourier transform beyond the fact that the Fourier 



16 



Table 2: Fractal dimension of fractals rotated by 45 degrees. 



Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdorff 


Dragon 


1.4339 


1.4535 


1.3839 


1.5760 


1.5236 


Fibonacci 


1.2525 


1.3479 


1.1742 


1.2499 


1.2465 


Gosper 


1.1389 


1.1240 


1.0627 


1.1213 


1.1292 


Julial 


1.1140 


1.0566 


0.8218 


1.1485 


1.0812 


Julia2 


1.1997 


1.1992 


0.9498 


1.2509 


1.2683 


Julia 3 


1.0959 


1.0189 


0.9748 


1.0076 


1.0043 


Koch 


1.2176 


1.2526 


1.1712 


1.2722 


1.2619 


Rabbit 


1.2783 


1.3657 


1.1106 


1.3586 


1.3934 


Tcrdragon 


1.2275 


1.0960 


0.9882 


1.2313 


1.2619 


Vicsck 


1.5059 


1.4368 


1.0311 


1.5436 


1.4650 



Table 3: Fractal dimension of fractals rotated by 90 degrees. 



Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdor 


Dragon 


1.4832 


1.4524 


1.4855 


1.5760 


1.5236 


Fibonacci 


1.3029 


1.3453 


1.2710 


1.2499 


1.2465 


Gosper 


1.1428 


1.1217 


1.1224 


1.1213 


1.1292 


Julial 


1.0940 


1.0538 


1.1612 


1.1485 


1.0812 


Julia2 


1.1468 


1.1972 


1.1226 


1.2509 


1.2683 


Julia 3 


1.0913 


1.0166 


1.2325 


1.0076 


1.0043 


Koch 


1.2487 


1.2497 


1.2141 


1.2722 


1.2619 


Rabbit 


1.3538 


1.3643 


1.2573 


1.3586 


1.3934 


Terdragon 


1.2410 


1.0949 


1.1132 


1.2313 


1.2619 


Vicsek 


1.4974 


1.4356 


1.3512 


1.5436 


1.4650 



Table 4: Fractal dimension of fractals scaled by 0.5 factor, that is, with 
halved size either in horizontal and vertical direction. 



Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdorff 


Dragon 


1.3439 


1.5337 


1.2109 


1.5760 


1.5236 


Fibonacci 


1.2267 


1.2895 


1.1491 


1.2499 


1.2465 


Gosper 


1.0669 


1.1078 


0.9326 


1.1213 


1.1292 


Julial 


1.0305 


1.0445 


1.0192 


1.1485 


1.0812 


Julia2 


1.0946 


1.1760 


1.0957 


1.2509 


1.2683 


Julia 3 


0.9157 


1.0133 


1.0706 


1.0076 


1.0043 


Koch 


1.1760 


1.2279 


1.0777 


1.2722 


1.2619 


Rabbit 


1.2057 


1.3776 


1.0768 


1.3586 


1.3934 


Terdragon 


1.1455 


1.2571 


0.8698 


1.2313 


1.2619 


Vicsek 


1.3349 


1.4803 


1.1157 


1.5436 


1.4650 



Table 5: 


Fractal dimension of fractals s 


icaled by a 


2 factor. 


Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdorff 


Dragon 


1.4948 


1.2672 


1.5031 


1.5760 


1.5236 


Fibonacci 


1.3586 


1.3354 


1.3830 


1.2499 


1.2465 


Gosper 


1.1608 


1.1318 


1.1538 


1.1213 


1.1292 


Julial 


1.1218 


1.0711 


1.0345 


1.1485 


1.0812 


Julia2 


1.1701 


1.2126 


1.1882 


1.2509 


1.2683 


Julia 3 


1.0017 


0.9677 


1.0492 


1.0076 


1.0043 


Koch 


1.2880 


1.2542 


1.0354 


1.2722 


1.2619 


Rabbit 


1.3042 


1.2495 


1.2831 


1.3586 


1.3934 


Terdragon 


1.2340 


0.8577 


1.3865 


1.2313 


1.2619 


Vicsek 


1.4046 


1.3318 


1.5783 


1.5436 


1.4650 
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Table 6: Fractal dimension of fractals translated by 10 pixels from the orig- 
inal position either in horizontal and vertical direction. 



Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdorff 


Dragon 


1.4396 


1.4524 


1.4137 


1.5760 


1.5236 


Fibonacci 


1.3045 


1.3453 


1.2832 


1.2499 


1.2465 


Gosper 


1.1260 


1.1217 


1.1029 


1.1213 


1.1292 


Julial 


1.0861 


1.0538 


1.0868 


1.1485 


1.0812 


Julia2 


1.1446 


1.1972 


1.1635 


1.2509 


1.2683 


Julia 3 


0.9820 


1.0166 


1.0952 


1.0076 


1.0043 


Koch 


1.2477 


1.2497 


1.0556 


1.2722 


1.2619 


Rabbit 


1.2723 


1.3643 


1.2078 


1.3586 


1.3934 


Terdragon 


1.1889 


1.0949 


1.0732 


1.2313 


1.2619 


Vicsek 


1.3767 


1.4356 


1.3319 


1.5436 


1.4650 



Table 7: Fractal dimension of fractals translated by 20 pixels. 



Fractal 


Box-counting 


Minkowski 


Fourier 


Proposed 


Hausdorff 


Dragon 


1.4412 


1.4524 


1.4267 


1.5760 


1.5236 


Fibonacci 


1.3022 


1.3453 


1.2843 


1.2499 


1.2465 


Gosper 


1.1297 


1.1217 


1.1158 


1.1213 


1.1292 


Julial 


1.0849 


1.0538 


1.0656 


1.1485 


1.0812 


Julia2 


1.1457 


1.1972 


1.0821 


1.2509 


1.2683 


Julia 3 


1.0025 


1.0166 


1.1686 


1.0076 


1.0043 


Koch 


1.2425 


1.2497 


1.1116 


1.2722 


1.2619 


Rabbit 


1.2722 


1.3643 


1.2115 


1.3586 


1.3934 


Terdragon 


1.1882 


1.0949 


1.0503 


1.2313 


1.2619 


Vicsek 


1.3688 


1.4356 


1.3679 


1.5436 


1.4650 



Table 8: Average error in the calculus of each fractal dimension calculus 
method in the fractals affected by classical geometrical transforms. 



Method 


Root Mean Square Error 


Box-counting 


0.2340 


Bouligand-Minkowski 


0.2249 


Classical Fourier 


0.3610 


Proposed method 


0.1269 



Table 9: Calculus of fractal dimension in fractals affected by 


a puntual noise. 




Fractal 


Box-counting 


Bouligand-Minkowski 


Classical Fourier 


Proposed method 


Waited value 


Dragon 


1.4461 


1.4623 


1.3712 


1.5404 


1.5236 


Fibonacci 


1.3070 


1.3512 


1.2580 


1.2494 


1.2465 


Gosper 


1.1408 


1.1357 


1.1737 


1.1275 


1.1292 


Julia 1 


1.0965 


1.0591 


1.1375 


1.1560 


1.0812 


Julia 2 


1.1529 


1.1997 


1.2146 


1.2395 


1.2683 


Julia 3 


1.0494 


1.0315 


1.1498 


1.1924 


1.0043 


Koch 


1.2503 


1.2542 


1.1289 


1.2384 


1.2619 


Rabbit 


1.2737 


1.3642 


1.2359 


1.3598 


1.3934 


Terdragon 


1.1871 


1.0949 


1.1162 


1.2313 


1.2619 


Vicsek 


1.3767 


1.4358 


1.3302 


1.5351 


1.465 
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Table 10: Average error in the calculus of each fractal dimension method 
applied to fractals affected by noise. 



Method 


Root Mean Square Error 


Box-counting 


0.2307 


Bouligand-Minkowski 


0.2244 


Classical Fourier 


0.3667 


Proposed method 


0.2228 



transform is an algebraic operation involving the manipulation of real val- 
ues, which turns it unsusceptible to rounding errors as occurs with the other 
estimation methods, mainly that based on geometric analysis. 

The Table [8] still shows that the proposed technique also is precise in 
the calculus of the fractal dimension, exhibiting the dimension value more 
similar to the theoretical value calculated by Hausdorff. 

Finally, we observe that the tables [9] and 10 demonstrate the robust- 
ness to noise presented by the proposed method. This is explained by the 
fact that the Fourier transform concentrates the noises into the last coeffi- 
cients, without a significant contribution to the dimension calculus, unlike 
geometrical methods which propagate the noise strongly. 



Table 11: Average computational time used by each method analyzed. 



Method 


Computational Time (s) 


Box-counting 


3.0394 


Bouligand-Minkowski 


4.2518 


Classical Fourier 


4.2482 


Proposed method 


0.4318 



Lastly, we verified the computational cost dispensed by each analyzed 
method, that, what the average time needed to calculate the fractal dimen- 
sion of a fractal represented in a digital image by each estimation method. 



The result is shown in the Table 11 The proposed method demonstrated 



a large advantage over the other methods being 7 times faster than Box- 
counting technique, the second faster algorithm. 



7 Conclusion 

This work proposed a novel method for the estimation of the fractal dimen- 
sion of two-dimensional fractal (objects or curves) contours. The method 
proposed is based on the exponential relation between the power spectrum of 
the Fourier transform of a signal and the frequency variable. It is a method 
that can be easily implemented in a computer and presents a significant 
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precision in the fractal dimension calculus as demonstrated by the results 
here exposed. 

The results showed that the proposed technique presents a precision in 
the calculus of the dimension greater than traditional dimension estimation 
methods like Bouligand-Minkowski, box-counting and classical Fourier di- 
mension. Besides, it was demonstrated the robustness of the method relative 
to the most common geometrical operation of rotation, translation and scale 
as relative to the presence of punctual noises. 

The precision of the novel method is explained by the inherent precision 
of the Fourier transform and the robustness of the technique is interesting 
once it allows the estimation of the fractal dimension independently from 
the geometrical position of the object and even when the fractal is affected 
by noises, which is very common in the practice, either due to noise added 
in the process of capturing or digitalizing of the object. In this way, these 
good results suggest the use of the proposed method in applications which 
yield to the calculus of fractal dimension in varied fields of the science. 
For a practical approach, the method could be used to estimate the fractal 
dimension of signals or temporal series (in the real or complex domain), 
founding a wide range of applications in science and engineering. 
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